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I. INTRODUCTION 



Since the discovery of high-T c copper-oxide superconductors [jlj, investigations of superconductivity and magnetism 
of d-electron based oxides with perovskite structure have gained new interest. Iron based oxides are not supercon- 
ducting, but exhibit interesting electronic properties such as disproportionation as in CaFeO^ or metalic properties 
as in SrFeOa. In addition, these ferrites often shown very peculiar magnetic properties. 

Almost stoichiometric S^LaFeaOg-a, (8 ~ 0) can be prepared under air by long annealing at 200°C or using 
electrochemical oxidation in alkaline solution Q , || , Q . All reported data of structural and physical properties show 
very similar features. 

The low-temperature Mossbauer data for the perovskite Sr 2 LaFe308.94 studied in Q are consistent with a 2 : 1 ratio 
of the two types of Fe cation with electronic characteristics fairly close to Fe 3+ and Fe 5+ , respectively. The Neutron- 
Powder-Diffraction data at 5K, show clear evidence that the low temperature phase exhibits an antiferromagnetic 
ordering More precisely, there is antiferromagnetic ordering among subcells of ferromagnetically-ordered spins. 
At room temperature only is observed an average value of the charge, which demonstrates the existence of a fast 
electron transfer. 

So, upon heating the compound passes from a mixed-valence phase through a first-order transition to the param- 
agnetic average- valence phase. Experimental data || for the specific heat show two peaks, the highest of which seems 
to be related to a first-order phase transition. The smaller peak, on the other hand, disappears as one fills in more 
oxygen vacancies in the sample, as was shown by the data of ref. Q] for the sample S^LaFeaOg, where the specific 
heat only shows one (large) peak at about the same temperature as in || . 

The inverse magnetic susceptibility for Sr2LaFe30s.g4 shows Q an abrupt change at about 200K with a minimum 
indicative of a transition to antiferromagnetic ordering; and above this temperature the susceptibility is almost field- 
indcpcndcnt; while below there is a pronounced dependence on the magnetic field which appears to be an intrinsic 
property of the phase, and suggests the possibility of a (weak) ferromagnetism, probably due to spin mis-alignment. 
The same sharp jump in the magnetic susceptibility is reported in refs. || and Q. 

According to previous structural studies this compound exhibits an almost cubic perovskite structure, the 
rhombohedral distortion being very small and not detectable by X-Ray Diffraction analysis. The magnetic structure 
was determined from the Neutron-Diffraction data 0] and is shown in figure 1. 
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The coupling between iron cations is of a superexchange-type via oxygen anions. However, the usual approximation 
consists in treating the system as that of iron cations interacting via an effective exchange integral. More precisely, in 
our case the coupling between two Fe 3+ ions is antiferromagnetic with the exchange constant denoted henceforth by[] 
J a , while that between Fe 3+ and Fe 5+ is ferromagnetic with the constant J/, while there is no coupling between two 
cations Fe 5+ , as this would occur along the diagonal of the (slightly deformed) elementary cube, and thus would be 
negligible in comparison with the other two couplings (see figure 1). In our particular case, these exchange integrals 
are obtained by fitting the experimental data with a model of localized spins with two kinds of nearest-neighbor 
interactions. 

Accordingly, we consider a model that describes a system of 6 localized spins, four of | and two of |, based on a 
Dirac-Heisenberg Hamiltonian with antiferromagnetic and ferromagnetic couplings J a , Jf, respectively and taking into 
account the corresponding biquadratic interactions, denoted by j a ,jf- The latter are known || to induce a first-order 
transition because of the strong change of the exchange interactions as a function of the interatomic distances, and 
also the associated change in the elastic energy of the material. The characteristic and very significant feature of 
such magnetic-transformation mechanism is the sharp change of the elastic constants of the crystal in the absence 
of any volume changes in the region of the antiferromagnetic transformation J5|. In our case, the exchange- induced 
distortion of the crystal is of an order of magnitude too small j|] , jj| to be responsible for the change of the order of the 
magnetic transition. In fact, it was pointed out in (^] and then shown in Q that the usual superexchange mechanism 
is fully capable of explaining the origin and magnitude of the biquadratic exchange interaction. Also Nagaev ( ||, 
section 2.7) studied the interplay between Dirac-Heisenberg and different kinds of non-Dirac-Heisenberg interactions 
and the role of the latter in changing the order of the magnetic transition viz in MgO, MnO and NiO. Moreover, since 
the spins in our case are much larger than i, more excited spin states start to participate in increasing the entropy 
of the system, and the two-level-system approximation no longer holds as the temperature increases. Therefore, in 
order to account for this effect, we must add higher-order terms to the Dirac-Heisenberg Hamiltonian. Accordingly, 
we shall add nearest-neighbor biquadratic interactions which, from the physical point of view could account for 
the effect of charge disproportionation occurring in our compound. Next, we show that indeed the biquadratic 
contributions to the exchange energy account for the first-order magnetic transition occurring in S^LaFesOg, and 
which is reflected by a sharp jump in the magnetic susceptibility and large peak in the specific heat. Consistently, we 
also predict an anomalous behavior for the sublattice magnetizations, i.e. the jump down to zero in the vicinity of 
the Neel temperature. Indeed, we obtain a good fit of the experimental data on the specific heat and the magnetic 
susceptibility in both the ordered and disordered phases for the exchange interactions jr- — 26K, it- ~ 6.5K ; and 
the corresponding biquadratic exchange integrals j a ~ 0.11 x J a and jf ~ 0.08 x Jf. 

Note that the antiferromagnetic exchange integral J a between two iron ions Fe 3+ in perovskite compounds such as 
LaFeCkwas also estimated by Anderson who found J a ~ 26K , and a close value, 2AK, was also found by Grenier 
et al. in CaFe2 0s with similar perovskite structure. 

Next, we compute the magnon spectrum within the approach of linear spin-wave theory. The spectrum contains a 
gapless antiferromagnetic mode and two gapped ferromagnetic branches. Unfortunately, to the best of our knowledge, 
there are so far no experimental data on the Inelastic-Neutron Scattering for the present powder sample. However, 
to account for the anomalous behavior of the sublattice magnetizations near the transition point, it turns out that 
non-linear spin-wave corrections are necessary. 

II. THEORY 

a. Hamiltonian The unit cell of our system is sketched in figure 1, and the Dirac-Heisenberg Hamiltonian for such 
system can be written as 

H DH = - Jf ^2 ^2 Sia ' S iP + Ja X! X! Sia ' S iP w 

where henceforth S,s denote the spins | of Fe 3+ , and | of Fe 5+ , respectively. denotes the sum over all pairs of 

{id) 

nearest- neighbor sites i,j; each Fe 5+ ion having six Fe 3+ ions as nearest- neighbors, while the nearest-neighbors of 



x The subscripts a and / refer to antiferromagnetic and ferromagntic orderings, respectively. 
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each Fe 3+ ion are three Fe 3+ and three Fe 5+ . The sum runs over the six different atoms in the unit cell, see 

a, 13 

figure f. 

Then to the Hamiltonian (Q) we add the following non-Dirac-Heisenberg Hamiltonian of nearest-neighbor bi- 
quadratic interactions 

h N dh = -if Y ( SiQ ■ s ^) 2 ~i a ^2Yl ( SiQ ■ s ^) 2 ( 2 ) 

a, p ot,P 

Note that we have adopted the convention that all exchange couplings J a , Jf,ja,jf are positive. 

Within the mean-field (MF) approximation it turns out that in fact the effect of adding these biquadratic interactions 
is to redefine the exchange integrals J and Jf, into temperature-dependent effective ones as follows 

J a — > Ja" = Ja+2ja(Sf (3) 

Jf — > Jf =Jf + 2j / (S) (s) 

where (S) and (s) are the spontaneous sublattice magnetizations on the sites Fe 3+ and Fe 5+ , respectively. The 
same is also true within spin-wave theory, see sect. 4. However within the linear-spin- wave approximation the effective 
couplings J^J ^ , are independent of temperature since in this case the spins S and s are respectively substituted 
for the sublattice magnetizations (S) and (s). 

As was demonstrated by Nagaev [||, all non-Dirac-Heisenberg interactions, biquadratic in our case, change the 
order of the magnetic transition from second to first, ft is readily seen in (^) that the effective exchange integrals 
decrease in magnitude with decreasing magnetizations because of a decrease in the order parameters (S) and (s). 
Inversely, a decrease in the effective exchange integrals in turn leads to a decrease in the magnetizations, i.e. there is 
a positive feedback effect. This brings about the change in the order of the transition as the ratio, j in the present 
case, between the non-Dirac-Heisenberg and the Dirac-Heisenberg interaction reaches a certain critical value. 

Henceforth, to avoid writing cumbersome formulae, we shall only give the expressions for physical quantities derived 
from the Hamiltonian (ij), but keeping in mind that all calculations are done using the effective exchange integrals 
Ja ff ,J} ff defined in @. 

b. Spontaneous magnetizations Within the MF approximation the order parameters (S 1 ), (s) are found to satisfy 
the following coupled self-consistent equations 

(S)=B S (zf3 [J a (S) + Jf(s)}) 

(4) 

(s)=B s (2zf3J f (S)) 
where f3 = -j^pf, and Bj(x) is the usual Brillouin function 

1 1 1 x 

Bj(x) = (j + -) coth(j + -)x - - coth - 

and z = 3, i.e. half the number of nearest neighbors of an atom. 

Upon making the substitutions (||) the self-consistent equations ([|) become more complicated and then they can be 
solved only numerically. The corresponding solution is shown in solid lines in figure 2, for the values of the exchange 
integrals obtained by fitting the experimental data on the specific heat and magnetic susceptibility, see sect. 3 below. 

c. Susceptibility We have computed the magnetic susceptibility both in the ordered and paramagnetic phases 
following the generalized mean- field approach of Smart 0. For this purpose, we use the fact that the spontaneous 
magnetizations given by eqs.(Q) exhibit the following antiferromagnetic-ferromagnetic arrangement at T < Tjy (see 
figure 1) 

(S 2 > = - <S 8 > = - <S B > = <S 6 > , (si) = -(s 4 ) (5) 
Thereby we obtain the following expressions for the transverse magnetic susceptibility, 

x Jf(2(S) + (s))> + 2J a (S)(s) 
2zJ f J a (S) 2 
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and parallel susceptibility^ 



2T- 



(2B' s (yS) + B'M) + 2z ■ (Ja + 4J/) B'MB' s (y$) 



X " aX T* + zJ a T.B' s (yl)-2z*jyB' s (y° j )B' s {yl) (7) 

where 

yl = (2zJ f (5)) f3, = z (J f (s) + J a (5)) (3 

and have introduced the conversion coefficient a = Ho(giJ-B) 2 N c /kB- With g = 2, and N c = ^,Na being the 
Avogadro number, as there are two molecules of the sample in the unit cell. 

It is worthwhile to note that the transverse susceptibility (|^) is a decreasing function of temperature. 

Now, as we are dealing with powder sample, the total magnetic susceptibility below the Neel temperature is given 

by 

X[T) = -X|| + -x± 

Next, in the paramagnetic phase the total magnetic susceptibility is found to be 

2 y + 2x + J a - AJf 

Xpm = a x — (8) 

z xy + J a x — ZJj 

where 

3k B T 1 3k B T 1 



x 



z s(s + l)' a z 5(5+1) 



One can check that as T — * Tjv, Bjx) — > J 'C?+ 1 ) ; so that eq.(Q) reduces to eq.(||), and at T = Tjy, we get = x\\ = 
Xpm ■ 

d. Entropy Within the foregoing approach the entropy (per atom) of the system reads 

= 2 [log (C a (/3A)) - (3 A (a)] + 4 [log (C s (f3B)) - (3B (5)] (9) 

where A = 2zJf (5), B = z [ J a (S) + J/ (a)], and Cj(x) — —^^rsj 2 > whose derivative is the Brillouin function 
given earlier. 

When plotted as a function of temperature, for the exchange couplings obtained below, the entropy (j^) increases 
with increasing temperature up to the Neel point where it exhibits an abrupt jump characteristic of a first-order 
transition, and then it saturates to a constant given by the configurational entropy in the paramagnetic phase. 

e. Specific heat The specific heat within the same approximation is given by^| 

^ = -2z- (J a ■ (5) d T (5) + J f ■ [(a) d T (5) + (5) d T (a)}) (10) 

for T < T N . 

Above T/v the mean-field approximation yields a zero specific heat since in this temperature range the order 
parameters and thereby the free energy vanish, see figure 4. 



2 Here the prime stands for the derivative of the Brillouin function _B(a;)with respect to x. 
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dr stands for the derivative with respect to temperature. 
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III. RESULTS AND DISCUSSION 



The Dirac-Heisenberg model including the biquadratic interactions studied here has been written for the ideal 
compound S^LaFesOg with no oxygen vacancies. On the other hand, the compound with the closest composition to 
the latter was studied experimentally by Wang et al. in M. Therefore, we believe that our theory is more suitable 
for fitting the experimental data of ref . [Q than those of or || , especially from the quantitative point of view. 
However, we also obtain good qualitative agreement with the authors of refs. ||, ||. 

A reasonable (numerical) fit to the experimental data Q on the magnetic susceptibility and specific heat given in 
figures 3 and 4, led to the following values of the exchange integrals and biquadratic interactions 

^ ~ 26K, i ~ 6.5K (11) 
k B k B 

j a ~ 0.11 x J a) j f ~ 0.08 X Jf 

As was mentioned in the introduction, the value of the antiferromagnetic exchange integral J a found here is in 
agreement with the ones obtained by Anderson jHj], i.e. J a ~ 26K in LaFeC>3, or by Grenier et al. |Tlj] , that is 2AK 
in CaFe205. 

For these couplings we find that the sublattice magnetizations exhibit a sharp drop at the Neel temperature 
indicating that the transition is of first order, see figure 2. This behavior is also reflected in the magnetic susceptibility, 
which exhibits a sharp jump at the transition, in agreement with the susceptibility measured by Wang et al. Q and 
also with the one reported by Zhou et al. ||, see figure 3. We see that the susceptibility curve obtained within 
the mean-field theory fits both qualitatively and quantitatively to the one measured by Wang et al. however the 
authors of ref. Q obtain a larger jump at the critical temperature. 

In addition, the jump in the entropy (|^) is indicative of a first-order transition, and yields the latent heat released 
by the system at the transition. Accordingly, for the exchange integrals found above, we find that the contribution 
of magnetic excitations to the latent heat is given by the product of the jump in the corresponding entropy at the 
critical temperature and the latter, i.e. Q(Tjv) = TV x AS ~ 2.9 KJ/mole. The experimental enthalpy found by 
Zhou et al. || is AH ~ 3.7 KJ/mole, which is the total enthalpy of the system. 



Consequently, the specific heat obtained within the mean-field approximation (10) diverges at the Neel temperature. 
In the magnetically ordered phase and around the transition we obtain a good fit of this to the experimental curve 
obtained by Wang et al. Q], see figure 4. As is well known the mean-field theory cannot be correct above the Neel 
temperature, for it predicts the absence of short-range order. In particular, in our case the specific heat computed 
within mean-field approximation drops to zero for T > Tm, since then the order parameters vanish, i.e. (S) = 
(s) = 0. Alternatively we have computed the contribution to the specific heat ( fio|) in the paramagnetic phase using 
the appraoch of high-temperature expansion, but this does not yield a significant contribution. We have also taken 
into account that we are only dealing with the contribution of magnetic excitations to the transition while leaving 
out the lattice component. 

To derive the Curie- Weiss law for the magnetic susceptibility in the paramagnetic phase, we note that the denomi- 
nator in (|8|) is quadratic in temperature, and thus leads to a hyperbolic function of temperature, as in the case of a 
ferrimagnet. However the high temperature asymptote to the hyperbola does have the Curie- Weiss law form, that is 

* = r%I < 12 > 

with 

C a 12, 9 a ~ -228K 
in fair agreement with the experimental result |2| , || 

C = 11.4, 9 = -250K 

In the paramagnetic phase, we can also estimate the average magnetic moment corresponding to the average valence 
of iron. Indeed, the relationship between the magnetic moment per atom and susceptibility 



yields for iron 



= y/8 x X (T - 6a) = 2.83 x ^ 



H(Fe) = 2.83 x ^ Ca ^ =s 5.66^ s 
which agrees with the value found by Battle et al. (|], that is 5.55/iB. 
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IV. SPIN WAVE THEORY 



/. Spectrum Within the linear spin- wave theory based on the Holstein-Primakoff representation 12 1 of spin opera- 
tors (-g expansion), we find three doubly-degenerate magnon branches, as shown in figure 5. The explicit expressions 
of the corresponding magnon energies are rather messy and we omit writing them here. Nonetheless, around the point 
r, k = (0, 0, 0), located at the center of the Brillouin zone we can write, 



V /p^-(fc2+fc2) +p l. fc 2 



hu 2 (k) = 

hj 3 (k) = Ply-{kl+kl)- Pl-kl 



(13) 



where the gaps Ai and A2 are given by 



Ai = ^ [2J a S S + J f (2S + s) 2 ] 



A, = 



^ (2J aS S + J f s 2 ) 



and the spin stiffness coefficients p\ y -,p\, etc., are (cumbersome) functions of the exchange integrals. Recall that the 
exchange integrals J a and Jf must be redefined using eq.(^) and taking into account the fact that in the lincar-spin- 
wave approximation the sublattice magnetizations (S) and (s) are replaced by their nominative values S = | and 
s = I, respectively. 

It is seen in figure 5 that the first of the spectrum branches, the lowest curve, is gapless and of antiferromagnetic 
type. Whereas, the second and third branches, the upper ones, represent gapped ferromagnetic modes. We have 
assumed here that the anisotropy is too small to produce a gap in the magnon spetrum at temperatures the latter is 
obtained. 

Therefore, the spectrum above shows that we have antiferromagnetic ordering at low temperature, together with 
a "weak ferromagnetic ordering" that starts propagating upon heating. As discussed in the introduction, it should 
be interesting to compare our results for the magnon spectrum with the experimental data, were Inelastic-Neutron 
Scattering measurements possible on the compound S^LaFesOt). This would also allow us to compare the spin 
stiffiness coefficients we have found here with the experimental ones, and then determine the effective magnetic 
moments of Fe 3+ and Fe 5+ at zero temperature, and thereby estimate their reduction by quantum fluctuations. 

g. Brillouin zone The Brillouin zone is a 3D— hexagon (see figure 5) whose edges are defined by 



— 7T , 7T 



-2tt 2tt 



-in 



< K < 



4tt 

T 



\ky\ 

V3 



h. Thermodynamic quantities The ground-state energy per site is given by 



Eg 



= H 



1 ^L8j2 



N c 2 



i/=i 



B.Z. 



d 3 k 



huj v {k) 



(14) 



where 



Ho = -zJ a ■ S(S + 1) - zJf ■ [S(s + 1) + s(S + 1)] ~ -0.16 eV 



and jf- 



-0.12 eV. 

Next, within linear spin-wave theory we have studied the temperature dependence of the internal energy, specific 
heat, and sublattice magnetizations. We have found that at low temperatures, the internal energy behaves as T , and 
that the specific heat behaves as T 3 , which is consistent with an antiferromagnetic ordering at low temperatures. 

On the other hand, as discussed in the introduction, the linear spin- wave theory turns out to be a poor approximation 
in the vicinity of the transition, as long as sublattice magnetizations are concerned. Indeed we find that (S) and (s) 
decrease linearly with temperature as this approaches the Neel point, and vanish at different temperatures, see figure 
2. There we also see that the linear-spin-wave theory yields, at zero temperature, quantum corrections to the magnetic 
moments of Fe 3+ ions but no corrections to those of Fe 5+ ions. Therefore, we think that further non-linear corrections 
should be taken into account so as to obtain the correct temperature dependence of the sublattice magnetizations 
especially as the temperature approaches the Neel point. 
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V. CONCLUSION 



We have modelled the first-order magnetic transition occurring in the perovskite S^LaFesOy using a Dirac- 
Heisenberg Hamiltonian including (nearest-neighbor) biquadratic interactions, whose origin could be related with 
the disproportionation of iron. We have been able, by fitting the experimental magnetic susceptibility and specific 
heat, to estimate the superexchange integrals for the antiferromagnetic and ferromagnetic interactions, as well as the 
biquadratic ones. The first of these integrals is consistent with the results of previous work on kindred compounds. 

Elastic-Neutron Scattering will be performed on this compound for checking the anomalous behavior of the sublattice 
magnetizations predicted here. Unfortunately, it is difficult to check up on the magnon spectrum we have computed. 
Nevertheless, the latter does confirm the magnetic structure determined from Neutron-Diffraction data by Battle et 
al. [| in the compound S^LaFezOs.^. 

Non-linear spin-wave calculations and Monte Carlo simulations are still under investigation. Finally, it should 
be very instructive to understand more in detail the effect of iron disproportionation, or more generally the charge 
redistribution, on the magnetic ordering. 
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Figure captions 



• Figure 1: The crystal and magnetic structure of S^LaFeaOg The large circles represent Fe 3+ and the 
medium ones represent Fe° + , with the arrows indicating the size and orientation of their spins; the smallest 
circles stand for intermediate oxygen atoms. The 2Sr 2+ : La 3+ cations which are disorderly located in the cubic 
centers have been omitted for clarity. The hexagonal unit cells of the two structures are commensurate, but 
the unit cell of the crystal structure is a triple one, while that of the magnetic structure is a primitive one, for 
the rhombohedral translations are absent in the magnetic structure. Each unit cell contains two formulae of 
S^LaFesOy, hence six iron ions, which form six magnetic sublattices below the ordering temperature. 



• Figure 2: The solid curves represent the temperature dependence of the sublattice magnetizations predicted by 
mean- field theory including biquadratic interactions. The upper (solid) curve represents the magnetization (S) 
of the Fe 3+ ions of spin |, and the lower curve is the magnetization (s) of Fe 5+ of spin §, obtained for the 
exchange integrals given in eq.([TT|) in the text. 

The dashed curves represent the corresponding sublattice magnetizations predicted by the linear-spin-wave 
theory, up to a temperature of 120 K . 



• Figure 3: The curve in balls represents the experimental magnetic susceptibility of ref. |4j], and the one in 
triangles is the suscpetibility measured by Zhou et al. The solid line is the magnetic susceptibility obtained 
from MF theory including biquadratic interactions. 



• Figure 4: The curve in balls is the experimental specific heat of ref. Q|, and the solid line is our theoretical 
result of MF theory including biquadratic interactions. 



• Figure 5: Plot of the magnon spectrum along the path ZYXMTY indicated in the Brillouin zone shown as 
inset. 
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